Direct observation of stress accumulation and relaxation in small bundles of 
superconducting vortices in tungsten thin-films 
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We study the behavior of bundles of superconducting vortices when increasing the magnetic field 
using scanning tunneling microscopy and spectroscopy (STM/S) at 100 mK. Pinning centers are 
given by features on the surface corrugation. We find strong net vortex motion in a bundle towards a 
well defined direction. We observe continuos changes of the vortex arrangements, and identify small 
displacements, which stress and deform the vortex bundle, separated by larger re- arrangements or 
avalanches, which release accumulated stress. 
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Among the most relevant topics in the physics of type 
II superconductors is the knowledge and control over flux 
line arrangements. Applications of these materials al- 
ways require fixing vortices by an adequate distribution 
of pinning centers. The characteristics and behavior of 
both, vortices and pinning centers, are critically depen- 
dent on the structural and superconducting properties of 
each sample. However, there are some important aspects 
which have general character. Usually, flux enters the 
sample from the edge in form of Abrikosov vortices, which 
are driven towards the interior by the Lorentz force due 
to the Meissner shielding currents. As first proposed by 
Bean [1 , motion toward sample center is hindered by the 
pinning centers, which tend to pin any vortex that passes 
by. A nonequilibrium state, the critical state, is created 
where the vortex density is largest in the regions where 
magnetic flux enters the sample. At zero temperature, 
the Lorentz force is balanced by the pinning force. The 
metastable landscape of vortices can be altered in several 
ways, as e.g. external current circulation or changing the 
applied magnetic field [THS]. A big effort, especially since 
the discovery of superconductors with higher critical tem- 
peratures, has been invested to understand the establish- 
ment of the critical state, and the nature and universality 
of the mechanisms behind the organization of vortices in 
pinning landscapes [THTT] . In some particular situations, 
isolated vortices, and its motion, have been studied with 
great detail, at high temperatures and/or at very low 
magnetic fields p!2l - [2Q] . However, there is still no direct 
information, at the level of single vortices, about the way 
the vortex system responds to alterations of the criti- 
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cal conditions at very low temperatures. For instance, 
starting with the system in a critical state, an increase 
of the applied magnetic field brings it to an overcritical 
situation, and vortices relax to a new critical state. In 
order to progress on the understanding of such a process, 
we have followed by Scanning Tunneling Microscopy and 
Spectroscopy, STM/S, the way that groups of around 
30 vortices change their positions when an initial criti- 
cal state is modified by small increments of the applied 
magnetic field (5H<CH. Experiments were performed at 
100 mK to avoid thermal creep, in zero field cooled con- 
ditions and in an amorphous thin film. Vortex positions 
are determined by linear extended pinning centers due to 
small surface roughness, which can be previously imaged 
using the STM in topography mode. The direction of mo- 
tion is conditioned by the distribution of pinning centers. 
Vortex bundles with roughly hexagonal vortex arrange- 
ments are gradually distorted, until they relax through 
small vortex avalanches. 

In previous work we have shown that superconducting 
amorphous W-based films, fabricated using focused-ion- 
beam [21 , can be studied with STM/S without a partic- 
ular surface preparation, showing perfect BCS behavior 
at zero field and vortex images with many vortices in a 
large field and temperature range [22| [23] . Here we use 
the same experimental set-up. The sample is 200 nm 
thick, and has been grown by sweeping the focused ion 
beam onto a Au layer previously deposited by evapora- 
tion on a Si substrate. Vortex pinning occurs through fea- 
tures appearing at the surface, namely linear depressions 
of 5 nm height, similar to those discussed in previous 
work [22 } [23 j . Vortex lattice STS images are built from 
the normalized zero bias conductance changes as a func- 
tion of the position in maps of 64x64 points. Each image 
is taken at a fixed field some minutes after changing the 
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magnetic field. Images require between 8 to 30 minutes 
each and are repeated several times at a given magnetic 
field to check for eventual vortex motion (creep) as a 
function of time, which was actually never observed [24 . 
We have studied up to four different regions, at different 
magnetic fields (between 1 T and a few T), finding the 
behavior discussed in the following. In Fig. 1(a), we show 
the topography of the area on which we focus here, where 
we have identified three main linear surface corrugation 
features, two in the upper left and right corners and an- 
other one in the bottom left part, marked in the figure 
with red dotted lines (see also [25]). 



by pinning through the surface features. At the bottom 
panel of Fig. 1(b), we show the vortex image obtained 
after increasing the magnetic field in steps of 0.004 T 
up to 1.2 T[25 . Vortices of the bundle between pins 
(e.g. orange and magenta) tend to move to the bottom 
right part of the image, following a path which is roughly 
parallel to upper and lower linear surface depressions. 
On the other hand, vortices close to the linear surface 
pinning centers (e.g. blue and green) move along them. 
Note that motion occurs in smaller steps separated by 
large jumps, and that there is some correlation between 
jumps. 




FIG. 1: (a) Topography (STM image) of a typical region of the 
W deposited film where red dotted lines mark small linear 
features at the surface. The lateral size of the image is of 
260 nm, with a corrugation of 5 nm. Tunneling conductance 
is of 3/iS, with 2.5 mV bias voltage. In (b) we show the 
initial (upper panel) and final (bottom panel) images of a 
series composed of 33 vortex images taken, in the area with 
the topography shown in (a), at 100 mK when increasing 
the magnetic field. The final image of the series shows the 
trajectories of four vortices, highlighted by colored points, 
(c) The trajectories of 17 vortices are represented in the x- 
y plane of the topography image (a). These 17 vortices stay 
within the scanning window during the whole sequence. Other 
vortices, not highlighted, jump in and out of the scanning 
window during the sequence. Lines joining points in (b) and 
(c) are guides to the eye. See Ref.[25] for more details. 

In Fig. 1(b) we show superconducting vortices in the 
same region obtained at 100 mK, after zero field cooling 
and increasing H up to 1.072 T. In the upper panel of Fig. 
1(b) we show the starting vortex arrangement. A vortex 
bundle showing distorted hexagonal order is found at the 
region between pinning sites. At the top and bottom 
part of the image, vortex arrangements are determined 




Angle to X axis (degrees) (nm) 



FIG. 2: (a) Histogram of the angles formed by the vector 
joining the initial and final vortex positions and the x-axis 
of the image (see Fig.l). The inset shows schematically the 
overall vortex motion, towards an angle of around 12° with 
respect to the x-axis, and the maximum overall deviations 
found with respect to this direction, (b) Histogram of the 
size of the displacements of vortices as defined in the 
text. Most frequent displacements 5li are of a few nm. 

All trajectories are shown in Fig. 1(c). To highlight 
the overall direction of vortex motion, we have plotted in 
Fig. 2 (a) a histogram of the angle 9 with respect to the X 
axis of the vector joining initial and final positions in each 
vortex trajectory. There is a net flux motion flow towards 
a preferential direction around Oq = —12°, i.e. towards 
the lower right part of the image at a path roughly defined 
by the upper and lower pinning lines, and close to a high 
symmetry direction of the hexagonal vortex arrangement 
within the bundle between pins. Clearly, a magnetic field 
gradient is found close to the studied region. When in- 
creasing further the magnetic field, the gradient increases 
and vortices move towards the direction of the maximum 
magnetic field gradient, being attracted to the lower right 
part of the image as given by = —12°. In Fig. 2(b) we 
show a histogram of the size of all displacements after 
each magnetic field step 6£i. 6£i is, for each vortex path, 
given hy 6 £i = \{ri — r^zt)!, where rt is the vortex po- 
sition in magnetic field step i and r^-l, the vortex posi- 
tion in previous step. Clearly, most displacements Sii are 
small, and sometimes large jumps, associated with large 
vortex re-organizations, occur. 

It is interesting to discuss the accumulated length 5i = 

i=n 

^ 5ii of each vortex path as a function of the magnetic 

i=l 

field. This shows the mechanism mastering the behavior 
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FIG. 3: Accumulated length 61 = ^ Sk (see text) vs. mag- 

netic field of four vortex trajectories. The color code is the 
same as in Fig. 1(c). In the upper left inset, the Fourier trans- 
form of the relative displacement 6li as a function of the ac- 
cumulated length 51 for one vortex (the one with blue points 
in the figure) is shown. Arrow marks the peak found at the 
inverse of ao = 40nm. The lower bottom inset shows the 
histogram of the inverse of the position of this peak in the 
Fourier transforms of the accumulated lengths of many vor- 
tex trajectories. The histogram peaks at the inverse of the 
expected intervortex distance a ^ 43-46 nm. 



of our system on its way along successive critical states. 
The result (Fig. 3) is that there is a succession of several 
small size steps separated by larger jumps. 

The repeated succession can be discussed by analyz- 
ing the Fourier transform of the dependence of the step 
size S£i for each vortex path as a function of Si. In the 
upper left inset of Fig. 3 we show an example. There is 
a peak close to a^^ = (40nm)~^. ao is close to the ex- 
pected intervortex distance a = (4/3) -^/^(^o/H)"^/^ for 
a hexagonal vortex lattice (a varies from 46 nm to 43 
nm between 1.072 and 1.2 T). This peak shows that 
the appearance of larger steps in 5ii occurs at regular 
intervals of the accumulated length Si of vortex paths. 
The peak is found close to a in many vortex paths, as 
shown by the histogram of its position (lower right in- 
set of Fig. 3). Thus, vortex paths in our experiment are 
modulated with lattice periodicity. This is remarkable, 
as translational order disappears within a few intervor- 
tex distances. Vortex motion showing some periodicity 
has been observed in vortex lattices driven at high ve- 
locities and/or high temperatures under different types 
of pinning potentials [H HSl [H HOl [H [27^ . Here, how- 
ever, the vortex lattice is static. The landscape of linear 
pinning centers in our samples allows the formation of 
nearly hexagonal vortex bundles in between pinning cen- 
ters at some magnetic field values. The orientation and 
position of these nearly hexagonal vortex bundles mimics 
the interplay between elastic energy, tending to organize 



the vortices in hexagons, and pinning energy, tending to 
fit them to the linear features at the topography. Such 
hexagonal very short range order possibly creates the ob- 
served regularity in the overall accumulated vortex mo- 
tion. 

Further microscopic insight can be obtained by ana- 
lyzing in detail vortex positions at each magnetic field 
step. In Fig. 4 we highlight a particular interval of six 
steps, between 1.112 T and 1.132 T. At the beginning of 
the sequence, the central vortex bundle shows a nearly 
hexagonal arrangement. Then, during the first four small 
steps (from 1.112 T to 1.128 T), vortices travel small 
distances and the central bundle becomes gradually dis- 
torted. Vortices with five or seven nearest neighbors ap- 
pear in the images, as shown in Fig. 4(b) (at 1.112 T, 
1.12 T and 1.128 T). In the next step, a dramatic rear- 
rangement occurs. It involves about ten vortices, which 
jump at the same time over larger distances Sii than in 
previous magnetic field steps (see Fig.4(a)). These collec- 
tive large jumps, small vortex avalanches, are the mecha- 
nism which governs the passage between two nearby well- 
differentiated vortex arrangements within a given pinning 
landscape. To our knowledge, these data are the first 
which evidence vortex avalanches at the scale of individ- 
ual vortices at magnetic fields as high as a Tesla and very 
low temperatures. Thus, the response of our system to 
the quasistatic increase of the external magnetic field is 
composed of two different processes, distorting quasielas- 
tic small displacements, and vortex avalanches, by which 
all or part of the accumulated stress is relieved. 

We can make a simple estimation of the elastic proper- 
ties of the bundles involved in our experiment. The pin- 
ning force Fp can be obtained from the magnetic field, 
assuming that Fp compensates the Lorentz driving force 
produced by the screening currents J which maintain the 
local magnetic field gradient [THS] . J = Jc, where Jc is the 
critical current density, which we estimate from macro- 
scopic measurements to be about 10^ A/cm^ here [21} [28]. 
and thus Fp ^ JcB. We obtain pinning forces of around 
1.1 10^ N/m^. The slope SFp/Sii gives, within a sim- 
ple Hooke's law approximation, an elastic constant. For 
example, we find a vortex lattice compression modulus 
Cii of around 300 N/m^, (taking Cn =FpRcdio/Sii with 
Rc being the radial correlation length [5 , 6 ) when we ex- 
amine the parts of the trajectories related to the regular 
vortex motion previously discussed (Fig. 3). The corre- 
sponding potential well is found to be around 5 K. Po- 
tential wells up to several tens of K are found, within this 
basic analysis, for motion involving avalanches of many 
vortices. 

Let us note that detailed simulations of superconduct- 
ing vortices within a pinning landscape given by point like 
pinning centers at zero temperature show a rather rich 
behavior, governed by vortex avalanches [29l[30]. When 
the density of pinning sites is low, easy channels for vor- 
tex flow appear, related to interstitial vortices moving 
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FIG. 4: In (a) we show the accumulated length 61 (see text and 
Fig. 3) of a set of vortex trajectories in a short magnetic field 
interval. Color code for identifying each vortex trajectory is 
the same as in Figs.l and 3. Dashed lines mark the magnetic 
fields where we have taken the images shown in (b). There, 
we show vortex arrangements at 1.112 T, 1.120 T, 1.128 T 
and 1.132 T. Vortices in the central bundle are marked by 
points and have been Delaunay triangulated. Green, blue and 
orange points show vortices with, respectively, five, six and 
seven nearest neighbors. Red dashed lines mark the pinning 
centers, as in Fig. 1(a). 



around their pinned neighbors, and at high pinning den- 
sities, vortices jump from pin to pin. In our experiment 
pinning centers are not point like but linear. Remarkably, 
vortices are able to stick to the pinning centers even if 
the pin is oriented perpendicular to the overall flow (as 
shown by the blue vortex of Fig.l during the first part 
of the sequence). This occurs at the expense of allow- 
ing jump and vortex re-orientations within the central 
bundle. 

In summary, we have studied vortex motion in the crit- 
ical state of a superconductor with linear pinning centers 
at high magnetic fields and in the absence of thermal 
excitation. Our experiment gives answers some relevant 
questions. Perhaps the most important is about the very 



nature of "avalanching objects" and has been expressed 
in the following way [7 : are they individual vortices, or 
flux bundles? are they rigid entities? Our experiments 
show that these objects are bundles of vortices, and that 
they are deformable entities. Even more, this deforma- 
bility is a fundamental aspect of the way bundles move 
between pinning centers. 
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EPAPS Figure 1. STM image of the area discussed in the paper. The black hne 
crossing a surface hnear feature (marked by red dotted hnes) is the profile plotted in 
the right panel. It shows the depression found at a linear surface feature. 
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EPAPS Figure 2. Vortex positions in the whole sequence of 33 vortex images. Red lines mark the position 
of surface linear features observed in the topography image of the studied area. In the initial frame, each 
vortex is marked by a colored number (labels for vortices in the same row have the same color). In the 
following frames, vortex positions are identified correlating consecutive images, starting with the vortices 
that move less, then their neighbors, and finally those with larger changes in position. Note that vortices 
appearing in the scanning windows are identified with a new color (e.g. yellow for vortices entering the 
image by the lower left corner in frame 11). 



